ECoG observations of power-law scaling in the human cortex 
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We report the results of our search for power-law electrical signals in the human brain, using 
subdural electrocorticographic recordings from the surface of the cortex. The power spectral density 
(PSD) of these signals has the power-law form P(f) ~ f~ x from 80 to 500 Hz. This scaling index 
X = 4.0 ±0.1 is universal, across subjects, area in the cortex, and local neural activity levels. The 
shape of the PSD does not change with local cortex activity, only its amplitude increases. We 
observe a knee in the spectra at fo c± 70 Hz, implying the existence of a characteristic time scale 
r — (27r/o) _1 — 2 — 4 msec. For f < fo we find evidence for a power-law with \l — 2.0 =b 0.4. 



The human brain is arguably the most complex struc- 
ture known to mankind and on the verge of starting to 
grasp its own inner workings. How do our brains com- 
pute? How fast do they compute? How do they store 
information? How universal is all of the above? Ever 
since the first electroencephalography (EEG) recordings 
in 1924, the study of the electrical activity of the human 
brain has focused on its prominent low-frequency fea- 
tures, in particular the excitatory and inhibitory rhythms 
at specific frequencies, like the a (10 Hz) and (3 (20 Hz) 
rhythms [1]. Traditional EEG studies are limited to 
/ < 100 Hz. The fundamental processes at the indi- 
vidual neuron scale suggest a role of higher frequencies: 
the time of flight of a spike along an axon, the synaptic 
neuro-transmittor diffusion time, the integration time of 
the dendritic arbor. They are all near or sub 10 ms [2[. 
Synchronization and correlations associated with them 
are expected to exist at least up to 1 kHz. 

Electrocorticographic (ECoG) recordings from the 
subdural surface of the cortex have recently made it pos- 
sible to examine the electrical activity of the human neo- 
cortex with finer spatial and temporal fidelity than EEG. 
An array of electrodes is placed directly on the surface of 
the cortex, see Fig. [TJ The absence of the skull and sur- 
rounding tissue increases the electrode voltage while the 
close proximity to the cortex means that ECoG records 
very local phenomena. For example, changes in the clas- 
sical ahj3 rhythms appear spatially uniform for a given 
set of tasks in EEG, but vary strongly spatially within 
the ECoG array for the same tasks [3[. 

The cortical surface potentials from sub-dural arrays 
reported in this study were obtained from 20 partici- 
pants receiving clinical monitoring for the localization 
of seizure foci prior to resection. Each was informed 
about, and consented to participate in, the University of 
Washington internal-review-board-approved experimen- 
tal protocol. The voltage was sampled at 10 kHz (2 sub- 
jects) or 1 kHz (18 subjects) using Synamps2 amplifiers 
(Compumedics-Neuroscan, San Antonio, TX) in parallel 
with long term monitoring (Xltek, Oakville, ON) from 32 
platinum electrodes, encased in silastic, in an 8x4 config- 
uration (4 mm in diameter, with 2.3 mm exposed, sepa- 
rated by 1 cm, center-to-center, Ad- Tech, Racine, WI). 

Our earlier studies H revealed the absence of distinct 



peaks in the power spectrum beyond / ^ 60 Hz. We 
hypothesized the existence of a power-law of the form 
P(f) — Af~ x at these higher frequencies, and named it 
the x-band/index, but the 1 kHz sampling rate truncated 
the signal at these high frequencies. The purpose of the 
study reported here was to determine, as accurately as 
possible, whether there is indeed such a power-law in 
the human cortical power spectrum, and how it might 
change with cortical activity (universality), by using a 
higher, 10 kHz sampling rate. 

Power-laws represent scale free behavior, the finding 
of which immediately evokes scale free networks, com- 
plexity, avalanches, and self-organized criticality (SOC). 
Unfortunately, many such networks and processes are 
not large enough or can not be monitored precisely long 
enough to establish the scale invar iance convincingly 
[BL m- The human brain is arguably the most com- 




FIG. 1: The electrode array locations are shown on a template 
brain for subject 1 (SI - purple, temporal) and subject 2 (S2 
- green, fronto-parietal) . Potentials of all 32 channels are 
measured simultaneously with respect to a scalp reference and 
ground before pairwise re-referencing. 
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Steps in Data Reduction 
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FIG. 2: Illustration of the steps taken to remove the amplifier 
roll-off and noise floors from the raw power spectra (green to 
dark blue to black) The amplifier built-in amplitude atten- 
uation was determined by scanning through 15-4000 Hz sine 
waves and fitting a smooth function to the attenuation at each 
frequency (light blue). A characteristic distribution of ampli- 
fier noise floors (current and potential noise) was generated 
by measuring the potential across an equivalent conformation 
of resistors (grey). The sharp line noise spikes at 60 Hz and 
its harmonics were excluded in our analysis. 



plex and largest network available and may provide the 
best opportunity to observe scale free behavior in a nat- 
ural setting. Each ECoG electrode measures the volt- 
age power spectral density from a specific cortical surface 
area associated with a specific set of functions. An elec- 
trode pair probes about 10 6 neurons, and each neuron 
has up to 10 4 synapses [2]. This has not gone unnoticed, 
and the neuroscience literature is awash with attempts to 
interpret experimental low-frequency data with scale free 
concepts and models* with only limited success and leav- 
ing many questions [7] . In this paper we firmly establish 
the existence of the power-law in the %-band and the ac- 
tual value of the scaling index, \. We obtain remarkable 
accuracy, particularly compared to many recent studies 
of power-law phenomena in nature [5[ . Our accuracy ap- 
proaches that of what was required and customary in 
equilibrium critical phenomena @. 

We characterized the power-law in the cortical spec- 
trum during a simple fixation task. The subjects fixated 
on an "x" on the hospital room wall 3 m away from the 
bed for 130 s (subject SI) or 190 s (subject S2), with 
their eyes open. The time-dependent voltage between 
each electrode and the reference electrode was measured 
at 10 kHz, digitized, and stored for spectral analysis. To 
further reduce the common mode noise from the envi- 
ronment, the digitized electrode voltages were converted 
digitally into a set of voltage differences between each 
near-neighbor pair of electrodes. This significantly re- 
moved the high-frequency common mode noise and make 
it possible to measure well above 100 Hz. For our 4 by 
8 array, the 32 individual electrode voltages produce 52 



near- neighbor voltages (Fig. [J). All of the measurements 
reported in this Letter are for near- neighbor voltages. 

We carefully characterized the amplifiers, their low 
pass filtering and their noise floor. These external fac- 
tors affect the power spectrum measurement dramati- 
cally. They mask and obscure the underlying power-law 
from the brain signal. Demonstration of the power-law 
would not have been feasible without these corrections. 
Fig. [2] shows these corrective steps in the data reduction. 
First, the power spectral density is calculated from the 
Fourier transform of the time- varying near-neighbor volt- 
ages in 1 s Hann- windowed epochs, overlapping by 0.5 s. 
These are averaged into uncalibrated spectra (marked 
green in Fig. [2]). They suggest a power-law shape, but 
the roll-off, apparent above 1 kHz, masks it. This roll- 
off does not represent a characteristic high-frequency in 
the brain. It originates in the amplifiers. We measured 
the gain versus frequency of the amplifiers independently. 
Their spectral bandwidth shape follows a classic low-pass 
filter (shown as blue in Fig. [2|). Removing the ampli- 
fier response from each uncalibrated spectrum leads to 
spectra (dark blue) that level-off at high frequencies, in- 
dicating that we hit a noise floor. This noise floor does 
not reside in the cortex either. We measured the am- 
plifier input noise for our amplifier-digitizer system inde- 
pendently. It is of order ~ 4 fiV per root Hz for every 
amplifier. Subtracting this from the spectrum leads to 
a power spectrum that remarkably tightly fits a straight 
line in the log-log plot (the black line see Fig. [2j) until 
at 500 Hz the signal disappears into the fluctuations of 
the amplifier noise floor. Future experiments with lower 
noise floor amplifiers will tell how high in frequency the 
power-law actually continues. 

Small variations and uncertainties in the amplifier 
noise floor C significantly affect the quantitative analysis 
of the power-law. For that reason we perform 3 param- 
eter fits, P(f) = Af-x + C. All such estimates for C 
are within the uncertainty of our independent amplifier 
noise floor measurements. 

Fig. [3] shows the combined spectrum, averaged over 
electrode pairs, for each subject. We exclude electrode 
pairs where one of the electrodes sits primarily on top 
of a vasculature, resulting in a much lower power in the 
signal and increased sensitivity to noise (4 channels for 
SI, 5 channels for S2). The inserts illustrate the quality of 
the power-law, the jitter around the average curve is more 
than one decade down from the signal, for all / < 200 Hz, 
beyond which the amplifier noise floor fluctuations start 
to kick-in. The exponent x an d the parameters A and 
C are estimated via a set of log-log type least-squares 
linear fits of the power spectral density between 80 and 
400 Hz, excluding harmonics of 60 Hz. This leads to x = 
4.03±0.1 for SI and \ = 3.96±0.1 for S2. The error bars 
are based on robustness against range shrinking as well 
as the deviations of the best fit with respect to the actual 
data across the entire frequency range 80 < / < 400 Hz. 

To test for universality, we also performed the exact 
same type of fits to each individual electrode pair spec- 
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trum. The histograms for SI and S2 in Fig. [3] over lap well 
and the variations within each can easily be attributed to 
stochastic external issues, such as variations in electrode- 
cortex distance and vasculature. A such, we reached 
the accuracy limit imposed by the experimental set-up. 
Within this accuracy, x 1S universal; it does not vary 
with subject nor specific brain areas. The mean value of 
the histogram exponents for SI is x = 4.06 (STD=0.10, 
N=48), and for S2 X = 3.94 (STD=0.13, N=47), consis- 
tent with the above fitting analysis on the two averaged 
spectra. We conclude that x = 4.0 ± 0.1 throughout the 
frequency range 80 < / < 400 Hz. 

Previous estimates of a power-law in the cortical spec- 
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FIG. 3: (A)-(B): Averaged electrode- pair power spectra 
(black) for subjects SI and S2. The red lines, the best power- 
law fits (see text) with \ = 4.03 (SI) and x=3.96 (S2), fit the 
data very well from 80-500 Hz. The spectrum before noise 
floor subtraction is shown in gray. The SI spectrum has a 
knee at at fo = 70 Hz and appears linear below fo. (C): his- 
togram of power- law fits from individual electrode-pair spec- 
tra, with mean values x — 4.06 (SI; green) and x — 3.94 (S2 
purple) . 



trum [8] focused on low frequencies. The averaged power 
spectrum of SI, in Fig. [3] shows a knee at fo ~ 70 Hz, and 
suggests a different power-law below fo. Both seem ab- 
sent in S2. The ah/3 rhythms are strongly pronounced 
in every channel pair of S2 (and are clearly visible in 
Fig. [3]). They obscure whatever power-law might be 
present underneath. Eight electrode pairs from SI lack 
aSz/3 rhythms. They may be absent in these local corti- 
cal areas, or so tightly synchronized that they cancel out 
in the electrode pair voltage difference. A simple minded 
P(f) — Af~ XL fit through their average from 15-80 Hz, 
yields xl = 2.57 ± 0.15 (N=8). However, such local fits 
are inherently dangerous. They are blind to the global 
properties of the spectrum. Our high frequency analysis 
of x is already a clear illustration of this. The amplifier 
roll-off and noise floors truncate our data only beyond 
/ ~ 400 Hz, but they affect the spectrum already at 
much lower frequencies (Fig. [2]). Similarly, the high fre- 
quency x = 4.0 power-law is already in-play below fo. 
Indeed, the global two-power laws form 

P(f) ~ A 



i+ 



= 4, shifts xl down to 



with constraint xl + Xh = X 
Xl = 1-99 (STD=0.21, N=8), with cutoff f L <C 15 Hz. 

The 8 SI channels form a dangerously small ensemble. 
We performed the same analysis on 1 kHz data for the 
exact same 2 minute fixation task from 16 subjects, with 
electrode arrays in the lateral frontal/temporal/parietal 
cortex. The spectra were determined and corrected for 
amplifier-digitizer frequency dependent attenuation in 
the same way as the 10 kHz data. We selected again 
only those channel pairs (N=116) that lacked the ahf3 
rhythms. Naive fits from 15-80 Hz (excluding 57-63 Hz 
line noise) of these corrected spectra yield again xl = 2.5 
(STD=0.4, N=116); while the same global two power- 
law form, shifts this down again to xl = 2.0 (STD=0.4, 
N=116) when setting = 1 and fo = 70. 

The uncertainty in the value of xl remains quite large. 
Even the existence of this low frequency power-law re- 
mains in question, until we collect more 10 kHz data. 
The analysis also requires a reliable decomposition of the 
low frequency ahf3 rhythms from the background. We 
demonstrated that principle component analysis (CPA) 
techniques can achieve this 0,0]. Fig. [4] illustrates this 
for 1 kHz data. Five subjects performed a finger move- 
ment task. A visual cue indicated when to move a given 
finger (repeatedly flex/extend one finger at a time; oppo- 
site side of body from grid placement), and the position 
of each finger was recorded using a dataglove (5dt; Irvine, 
CA). Samples of movement spectra of each type were cal- 
culated from the FFT of 1 s Hann- windowed epochs cen- 
tered at the maximum displacement of the finger during 
each flexion. Rest spectra were calculated from epochs 
in which there was no hand movement. The spectra were 
corrected for amplifier-digitizer frequency dependent at- 
tenuation. Differences between variations in ah/3 versus 
high frequency allows a PC A removal of the ahf3 peaks 
from all the power spectra. We were able to identify elec- 
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FIG. 4: The average shift in power spectral density in elec- 
trodes after the PCA decomposition in the hand cortical area 
during finger movement, for subjects 3-8; 5 electrode-pair 
channels each (30 total). (A) demonstrates that movement 
(red) increases the overall power by a factor 1.8 while preserv- 
ing the shape of the rest (blue) spectrum, f > fo power-law 
fits (grey) are consistent with \ — 4- (B) Electrode pair chan- 
nel locations (interpolated) across all subjects, projected to 
the left-hand side. (C) Remaining spectra after dividing out 
a Lorentzian, 1/ (l + (//70) 2 ) to illustrate consistency with 
the two-power laws form with \l — 2. 

trode pair channels which demonstrated significant shifts 
in power during specific activity (single finger movement 
type) versus rest. A subset of five such channels was cho- 
sen naively (based upon significance rather than being 
hand picked) from each subject, and the average spec- 
trum for the specific movement type causing change (the 
"active" state) was compared to the average rest spec- 
trum (the "inactive" state), after the ahf3 peaks had 
been decoupled and largely removed. 

Fig. [H addresses several issues: Earlier [3|, 0] , we ob- 
served changes in the x~band during behavior tasks and 
hypothesized how those can be used directly to quantify 
activity in the brain in a variety of practical settings; 
e.g., the total power in the x~ band increases with ac- 
tivity versus rest. Fig. [4] demonstrates that this increase 
in spectral power with activity extends over all frequen- 
cies. Moreover it strongly suggests that the shape of the 
spectrum is preserved (universality). The above xl ~Xh 



two-powerlaws form fits the shape well; with fo ^ 70 Hz, 
X — 4, and xl — 2, within 10%; a level of accuracy to 
be expected within the limitations of the 1 kHz nature 
of this data. The active/inactive power ratio R (between 
the amplitudes A) is unlikely a universal number. In this 
data set its geometric mean is equal to R = 1.81 with a 
variation of order 0.34 (N=25). 

What does all of this teach us about computations and 
correlations in the brain? EEG and ECoG voltages repre- 
sent the superposition of the electric current dipole fields 
generated by the very large collective of nearby neurons 
and their associated ionic channel currents, propagated 
through the complex mass of ionic extracellular liquid 
and neuronal and glial membranes. This is a quite com- 
plex phenomenon and not yet well understood quanti- 
tatively. Moreover, the underling neuron computational 
issues remain in flux, e.g, understanding the relative roles 
of computations at the level of dendritic trees versus 
those at larger length and time scales associated with 
the connectivity the network of neuron. 

If the value of the high-frequency exponent x — 4.0 ± 
0.1 were distinct from an integer, we could have reported 
safely to have observed SOC type complex scale free be- 
havior in the cortex. The uncertainty in xl — 2.0 ± 0.4 
still leaves room for this. Perhaps SOC behavior (if it 
exists) is only expressed in more subtle ways in ECoG. 
In the case xl = 2, the spectrum is well described as a 
product of two 1/(1 + (f / f c ) 2 ) Lorentzian shapes. These 
can originate without any SOC complexity: such as white 
noise with two filters; as the product of an exponential 
decaying correlation function and a temporal form factor; 
as two processes with definite correlation times; etc. 

Our results place definite constraints on future brain 
modeling. ECoG spectra definitely scale as P ^ f~ x 
across all 70 < / < 400 Hz. The value x = 4.0 ± 0.1 is 
universal, across subjects, areas in the cortex, and local 
neural activity levels. The knee in the spectra at fo ~ 
70 Hz, implies the existence of a characteristic time scale 
r = (27r/o) _1 — 2 — 4 msec, probably originating at the 
neuron length scale and surviving the coarse graining to 
the 2 mm electrode size length scale. 
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